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Q 

Three  wars  nf  wind  anil  temperature  data  drtrnniniil  from  limindr  experiments  al 
^six  nitis  an*  utilized  in  tin-  calculation  of  the  Itirlmidxon  ntnulH-r  JIM.  Tin*  results  show 
that,  for  a  critical  Iticlianlsoii  iiiiiuImt  of  I.  then-  is  i|uitc  often  a  turlnilcnt  region  around 
Ho  km  at  all  latitudes.  Furthermore.  tin-  winter  polar  mesospheric  data  indicate  a  more 
intcnsi-  turlndence  than  do  tlir  summer  data.  Then-  also  ap|M-ars  a  region  of  marginal 
staliility  ( 111  «•  I)  at  nr  slightly  shove  the  stratnpaiisc  which  is  more  often  pn-sent  during 
the  w  inter  season  than  in  the  summer.  I  altitudinal  and  seasonal  ilistrihntionsof  Itichardson 
niinds-r.  wind  shears  and  estimates  of  heating  rates  due  to  viscous  dissi|Hition  of  turhiih-nt 
kinetic  energy  are  given. 


I.  Theory  and  Analysis 

In  »  stably  stratified  atmosphere.  tlu-  neccssiiry  etinditioli  for  tile  litainte- 
nance  of  t orlinletlee  against  liiiovant-y  forces  Is  that  the  flux  Kiehnrilson  number 
In-  less  than  some  eritieal  niiiiilier  /i’f  •  where  lttt  is  the  eritient  flux 

Richardson  titiniln-r.  IT  this  flux  Biehnrdsott  iinmlM-r  eattnoi  In-  itieasttml  since 
it  rei|ttires  the  in  *itn  simultaneous  measurement  of  the  turbulent  transport 
coefficients  of  heat  ami  momentum  as  well  as  the  lociil  atmospheric  temperature 
and  wind  shear,  then  a  measure  or  stability  or  instability  may  In-  provided  by 
examination  of  the  ordinary  Itichardson  numlN-r  Hi  defined  ns 

/.*/  'j\r:  ?y.  ■  r\  T\?r  ?z\*  <d 

where  >j  is  the  neeeleration  of  gravity.  T  is  the  atmospheric  tcm|N-rnturc  (  K), 
/’is  tin-  adinlmtie  la|cs-  rate  (^lo  km'1).  X  is  the  vertical  coordinate,  and 
/’  is  the  amplitude  of  the  mean  horizontal  wind  velocity,  such  that  dl’j?Z 
is  the  vertical  wind  shear,  defined  by 

\?r  ?z?  \?r,  ?z\-  | ?rt?zp. 

and  where  r  and  7  are  the  horizontal  directions. 

The  amplitude  of  the  critical  Itichardson  numlN-r  Hir  necessary  for  the  onset 
of  turbulence  in  a  stably  stratified  atmosphere  has  In-eti  theoretically  shown 
to  In-  /•’/",  <  1  and  has  In-cii  ex|N*rimentnlly  verified.  WoihIs  |l|  has  shown 
t  hat  when  a  region  has  In-cit  placed  into  a  rood  it  ion  of  fully  devclo|N-d  t  urhulent 
flow,  it  may  Ik*  maintained  in  this  ttirhnh-nt  flow  under  the  condition  Hi  <,  1. 
and  since  for  this  study  we  do  not  have  a  time  and  space  history  of  the  tern- 
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peraturc  and  wind  fields,  we  shall  assume  that  the  observations  where  Hi  <£  1 
arc  indicative  oflocal  turbulence. 

Thus  given  this  condition  for  turbulence  (Si  1),  let  us  now  consider  the 
turbulent  energy  balance  equation,  which  may  be  expressed  as 

e/ATni  =  (l  —  «  Hi)  (dUfiZ)2  (2) 

where  e  is  the  rate  of  viscous  dissipation  of  turbulent  kinetic  energy,  Km  is 
the  coefficient  of  turbulent  momentum  transfer,  Kh  is  the  coefficient  of  t  urbulent 
heat  transfer,  and  a  =  Kh/Km,  assumed  unity  here. 

Utilizing  the  Heisenberg  [2]  relation  for  Km  and  assuming  an  “inertial" 
sub-range  ovc-  the  turbulent  energy  spectrum,  an  expression  for  e,  valid 
within  uncertainties  of  2,  may  be  derived  using  one-dimensional  theory, 

e  *  0.250  J2  [(1  -  Hi)  (S  £//0£)2]3/2  (3) 

where  l  is  the  large  scale  end  of  the  “inertial"  'spectrum.  To  limit -the  eddy 
wavelength  to  values  smaller  than  the  wavelengths  of  the  wind  shear,  we 
shall  use  a  scale  size  l  of  500  m  as  the  value  of  this  large  scale  limit.  This  will 
give  a  wavelength  of  3.14  km,  which  for  the  most  part  is  small  compared 
with  the  wind  system  vertical  wavelengths  of  5-20  km  usually  reported. 


RICHARDSON  NUMBER  (Ri) 


Fig.  1.  The  Churchill  (00°  N  but.)  winter  and  summer  Kielmrdson  numbers.  The  average 
values  are  shown  by  the  dashed  line,  and  one  immediate  sample,  ns  dated,  is  given  by  the 

solid  line. 
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Fixing  the  large  scale  limit  can  introduce  some  error  in  the  determination 
of  e,  since  by  theory  (Eq.  (3))  e  and  l  are  closely  coupled.  However,  this  method 
allows  A'),,  which  is  not  calculated  here,  to  vary  as  the  rate  of  dissipation  and 
thus,  within  the  limit  of  this  scale  size,  affords  a  more  self-consistent  deter¬ 
mination  of  £  than  would  be  accomplished  by  predetermining  A^ra  rational 
procedure  which  would  lead  to  serious  inconsistencies  due  to  the  aforementioned 
interrelationship  of  e,  A'h  and  l.  The  errors  generated  by  holding  l  constant 
would  be  to  overestimate  e  in  regions  of  high  shear,  and  thus  lower  Richardson 
numbers,  and  to  underestimate  e  in  the  region  of  Hi  near  unity,  or  low  shear  region. 
The  occurrence  average  of  e  and  Hi,  to  be  given  later,  should  be  considered 
conservatively  because  of  the  small  data  sample. 


2.  Results 

Murphy  et  al.  [3]  used  eight  years  of  wind  and  temperature  data  (1959-1967), 
as^determined-from  grenade  data-  [4], :  in  their  analysis,  bht  because  of  large 
errors,  particularly  at  the  uppermost  data  point  in  the  years  up  to  1965,  the 
results  of  tho  analysis  are  open  to. some  question.  We  shall  limit  this  analysis 
to  the  years  1965,  1966  and  1967;  even  so,  the  general  points  raised  in  the 


l’ig.  2.  The  Hichardson  numbers  for  the -Wallops  Island  (38°  A’)  data.  Tho  identification 
of  the  curves  is  tb>  same  nsin  Fig.  1. 
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eight-year  analysis  are  repeated  in  this  three-yeatanalysis.  The  winter  northern 
(60°  N)  mesosphere  shows  a  more  intense  turbulence  than  docs  the  summer 
northern  mesosphere  (Fig.  l)  ,snd  the  summer  data  are  somewhat  sporadic  in 
the  occurrence  of  turbulence  above  80  km,  while  the  winter  data  show  turbulence 
as  being  present  100%  of  the  time  at  some  altitude  above  75  km.  The  mid- 
latitude  data  (38"'  X)  (Fig.  2)  show  little  turbulence  for  the  years  1965  and 
1966,  but  more  frequent  occurrences  for  1967,  and,  as  an  overall  average, 
approximately  two  thirds  of  all  data  have  some  turbulence  above  80  km.  The 


Fig.  3.  The  Richardson  numbers  for  the  Xntnl  (f>°  S)  data.  The  identification  of  the 

curves  is  the  same  as  Fig.  1. 

equatorial  region  (5°  S)  shows  fairly  intense  turbulence,  present  most,  of  the 
time.  The  region,  at  or  slightly  above  the  strdtqpausc  (Figs.  1,  2  and  3)  shows 
sections  of  marginal  stability  or  even  turbulence  (lfi  ^  1,  and  occasionally 
<£  1/4),  and  ns  theoretically  demonstrated  by  Hines  [6],  these  regions  may 
correspond  to  a  critical  layer,  signified  by  Zfi  <£  1/4,  which  results  in  ampli¬ 
fication  of  the  background  wind  system.  Alternatively  if  the  condition  7 ’»'*’  <£  1 
signifies  a  turbulent  region,  ns  hypothesized  here,  the  radiation  of  short  internal 
waves  from  this  region  of  turbulence  arises  from  the  growth  of  wave  modes 
in  the  interactions  among  the  turbulent  components  as  noted  by  Phillips  [0J. 
’In  either  event,  it  is  possible  that,  these  stratospheric  regions  (Hi  <[  1)  are 
a  source  of  wave  generation  or  amplification  which  propagate  through  the 
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miBOKpheiv  to  ionospheric  heights  and  provide  the  stratosphere  -sonoapherie 
coupling  as  noted  by  ionospheric  pbysieietc. 

The  win^er-Kii’nmor  averages  of  the  rat“  of  viscous  dr  sipstion  for  the  three 
latitudes  are  given  in  Pigs.  4,  f>  ord  (i  as  a  function  of  aihtude.  Tin*  values  of 
r  ujHin  occasion  go  as  high  as  ~i<(*  ergs  g" 1  s  •  and  fairiv  often  up  hr  and  gt cater 
than  10s  ergs  g-1  s-1,  in  rt gions  of  unusually  h'gh  shears  (~10  *w  */.  parti-, 
eiilarly  at  the  northern  latitudes.  Because  of  tie*  p'ohhbie  eiror  genes  at  ion 
of  r  discussed  earlier,  values  of  ^  i  \  111®  an*  •  n.oved  from  the  cateulatmn 
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of  the  mean  hunting  rate ,  also  shown  in  Fig.  4.  6  ami  6.  Another  factor  to 
consider  in  removing  these  values  is  the  low  size  of  the  data  sample,  where 
a  single  value  off  >  10*  in  .t  sample  of  li  would  overwhelm  the  more  statist¬ 
ically  significant  sample  for  i  <*  10*  ergs  g~>  s  '1..  It  is  quite  apparent  that 
the  average  wint  -r  northern  mesosphere  has  a  rate  of  heating  which  is 
equivalent  to  or  higher  than  that  ealeulated  from  the  subsidence  and  m  cm., 
bination  of  atomic  oxygen,  and  that  of  Johnson  and  Gottlieb  [7 j  who  cab 
culated  a  heat  source  due  to  the  meridional  wind  distribution  which  was 
required  to  balance  the  lack  of  symmetry  in  the  solar  heat  input  for  solstice 
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conditions,  and  which  would  then  provide  the  ct.cigy  necessary  to  maintain 
the  warm  polar  mesosj  oere.  The  mid-latitude  heating  rate  is  obviously  much 
lotvei  than  that  of  the  northern  latitudes  ot  equatorial  region,  and  is  approxi¬ 
mately  of  the  same  amplitude  for  the  yearly  average  as  that  determined  by 
Zimmerman  and  Rosenlierg  [8J  from  wind  trail  analyi i3  at  approximately  95  km. 
The  equatorial  region  also  indicates  large  heating  rates,  similar  to  those  of  the 
northern  latitude.. 
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3.  Summary 

Thin  study  shown  the  atnionphere  above  80  km  in  usually  turbulent,  with 
the  occurrence  of  turbulence,  at  the  moment,  an  some  non-predietive  function 
of  latitude  and  season.  The  average  heating  of  the  atmosphere  by  the  turbulent 
flow  system  is  of  the  order  of  or  greater  than  the  heating  due  to  euv  and  Sehu- 
mann-Runge  energy  absorption,  and  apparently  com|>enaateH  for  the  lack  of 
symmetry  in  the  solar  heating  input  at  solstice  conditions. 

The  region  around  the  stratojwuise  shows  conditions  which  an1  sporadically 
favorable  for  wind  system  amplification:  however,  many  more  comparisons 
of  these  data  with  ionospheric  mesospheric  measurements  and  analysis  must 
lie  performed  to  ascertain  this  {tossihlc  influence. 
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